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We study numerically collisions between protons and H2 molecular ions at intermediate impact 
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sections are computed. The results are compared with: a) the standard one electron-two nucleon 
scattering, and b) the quantum mechanical treatment of the H + — H% scattering. 
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I. INTRODUCTION 

Ion-atom collisions represent one of the main fields of research in atomic physics, both experi- 
mental and theoretical. Currently, there is a great deal of studies about the collisions with electron 
transfer between ions and oriented atoms [i.e. with a preferred sense of circulation of the electron 
around the nucleus) and also between ions and aligned atoms (where the probability distribution of 
the electron is not spherical). The interest has been triggered by the development of experimental 
techniques which allow to prepare atoms in well defined states of low [jjj or high quantum numbers 
(n, I, m) H and, after scattering, to measure the final state of the system ||. The very narrow phase 
space volume sampled by the electron allows a detailed study of the physical mechanisms occurring 
during the impact, which would not be otherwise transparent due to the averaging over the entire 
space of electronic configurations. Consequences of first order effects-spatial overlap and velocity 
matching-have been extensively studied with the above mentioned techniques, interpreted in the 
light of classical mechanics jl|,^] and sometimes compared with exact quantum calculations ||[?]|. A 
recent brief review about the latest developments in this field is given by Schippers ||, and some 
even more recent experimental and theoretical works can be found in ^J. 

Less attention has been paid to processes which involve more than one target nucleus: some studies 
of collisions between H2 and bare ions are reported, for example, in |p||Io|-|l3t , and the scattering 
H + -Hf is the subject of [dj. 

Calculations on these systems using quantum mechanics is far from being straightforward. By con- 
trast, the application of classical methods, and in particular of the Classical Trajectory Monte Carlo 
method (CTMC) presents several advantages: The numerical complications introduced by solving 
the equations of motion for a few more particles are negligible. One may think, in comparison, to 
the problems which arise when attempting of solving the Schrodinger equation for three particles 
instead of two. 

Since the original work of Abrines and Percival ^| the CTMC method has been one of the most 
successful techniques for studying the scattering between heavy charged particles at intermediate 
impact energies. Starting from the simple H-H + processes and from the calculation of simple total 
cross sections, the method has been refined so that nowadays it gives many detailed informations: 
Differential cross sections Jl6|| , which are ordinarily measured in experiments, and state-to-state 
transitions Jit| , which are especially useful for the research and development on thermonuclear fu- 
sion. The CTMC method has also been applied to more complex systems: Collisions involving more 
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than an electron JlH , requiring the inter-electronic correlation, have been, and still are [^Jl9|. a 
great challenge. 

To our knowledge, until now CTMC has not been applied to polynuclear targets: In this work we 
do an investigation on this kind of process using the simplest target, the H2 molecular ion. We aim 
to understand how and to which degree the complex structure of the target modifies the results of 
the scattering with respect to the standard two nucleons-one electron case. 

In presence of a diatomic target the question of interference between the two scattering centres raises, 
and it has been faced within the quantum-mechanical formalism in Refs. |l(j,|ll],[L4]] . In classical 
calculations one cannot speak of interference as in quantum ones; however a certain modification of 
the results due to the presence of a second target is likely to be expected. 

The first problem one has to face is to obtain an equilibrium electron distribution function for 
the H2 molecular ion. In Section |n] it is described how this has been dealt with in this work. 

The results of the numerical simulations, expressed in the form of cross sections, are compared 
with the similar results for two nucleons scattering, and the differences are examined in Section III . 



II. THEORY 



In any CTMC calculation an important role is played by the choice of the initial conditions of 
the electron in the phase space. It is well known that, in the case of the ground state of hydrogen, 
extracting electron coordinates from a microcanonical distribution, f(r, p) = 6(E — Ei s )/8ty 3 , yields 
the correct quantum mechanical momentum distribution [19]: 

where _Ei s = —0.5 au (atomic units will be used unless otherwise stated), and N the normalization 
factor. Obviously this relations holds also for hydrogen-like ions. On the other hand, the radial 
distribution is not reproduced satisfactorily. While no classical method can reproduce at the same 
time both the exact momentum and spatial distribution, some methods have been devised which 
allow to exactly reproduce the latter distribution at the expenses of the former (so called CTMC-r 
method, opposed to the CTMC-p method-see |m| ), or to yield an approximate-but rather good- 
description of both distributions (for more about the subject see the references quoted in [n6| or 
ElJ). In this work we have approximated the quantum mechanical electron wave function by a 
Linear Combination of Atomic Orbitals: 

m = <(r)+<(r) (2) 
v2 

where A and B refer to the two protons which are placed initially with null velocity along the z axis, 
at a distance ±20 from the origin, with zo kept equal to 1 au, in accordance with the true equilibrium 
internuclear distance. With this choice the molecule has a cylindrical symmetry with respect to the 
z axis. From a Fourier transform of Eq. (2) we get the momentum distribution function: since 

4>t(r) + i>?s(r) = i>is(r-r )+ip la {r + r ) , (3) 

where ro = (0, 0, Zq) , we obtain 

HP) = j^yt J dr ^(r) e- lp - = -L (^^.(p) + e ip **>i s (p)) . (4) 
The probability density is better expressed in cylindrical coordinates 

g(p) = 2np r \i>{p)\ 2 = N' J r 2 cos 2 (p z z ) , (5) 
v r ~r Pr ~r Pz ) 

where Pz,p r are the projections of p along the z axis and the radial direction in the x-y plane. N 
is a normalization factor. 

The probability density q is similar to a hydrogenic distribution (see Eq. 1) but for the factor 
PrCOS (pzZo). This means that electron velocity is preferentially found within the x-y plane, where 
p z ^0. 

The couples p r , p z are picked up within a range {0,p m ax) and generated according to the dis- 
tribution of Eq. (5) with a rejection technique; p m ax is chosen great enough so that o(p r ,Pz) is 
negligible for pi + pi > Pmax- in the computations p ma x = 4.5. 
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The choice of the spatial position of the electron follows a similar route. The position r is 
constrained to satisfy 



_£ i i_ = E + 

2m e r — r |r + r | H i 

where ro = (0, 0, Zo) and E„+ ~ —1.1 au is the experimental quantity, r is characterized by the 
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three numbers: the radial distance r = |r|, the polar angle 9 and the azimuthal angle (f>. Again, r 
and cos(#) are chosen with a rejection technique within the ranges (0, r ma x) and ( — 1, 1) respectively. 
<f> does not explicitly appear in Eq. (6) and may be chosen uniformly in the interval (0, 2n). 
The chosen distribution is not expected to be a stationary one for the classical system, so it is of 
interest to see how much it varies with time. In Figure 1 we plotted the contour of the distributions 
Q{pz,p r ) g(z,r) at t = and t = 5, in absence of the projectile, from which one may see that the 
difference is not small but the essential features still remain so the choice of this distribution appears 
to be justified. 

The projectile initial parameters are the velocity v, the impact parameter b with respect of the 
centre of mass O of the molecule, the azimuthal angle i? with respect to the molecule axis, and the 
initial distance d from O. A sketch of the scattering configuration for the coplanar case (all the 
nuclei lying on the same plane) is shown in Figure 2. The projectile is then rotated out of the plane 
by a random angle between and 2n. 

After having initialized all the four particles, the corresponding equations of the motions are 
numerically integrated in time until the nucleons are well far apart. In the computations d is not 
a critical parameter as far as it is great enough to allow the target molecule and the projectile to 
be considered initially as non-interacting. After some trials a value of d ~ 20 au was found to be 
reasonable. With given values for b, v and a number TV of runs have been carried on, varying 
only the electron initial conditions. At the end of each collision one may find one of the following 
situations: 1) the original molecule remains intact; 2) it may be dissociated but the electron is 
still bound to one of the two nuclei; 3) the electron is bound to the projectile (we call this case 
"charge transfer"); 4) the molecule may be broken and the electron be ionized. Each process i (with 
i = l,2, 3, 4) happens Ni times over the total TV runs, with corresponding probabilities 

The standard deviation error of Pi is 



NN, 

Cross sections are computed by integrating over the impact parameter b 

Gi{v,$) = 2tt ( dbb Pi(b) . 
Jo 

where Pi (b) ~ for b > b ma x ■ Still, the cross section of Eq. (9) may be averaged over $: 

1 f 2n 
2n Jo 

This latter integral has been evaluated by a simple trapezoidal rule using the values -&k (k=l,...,n) 
for a finite set of angles. 



III. RESULTS 



The runs have been performed for 0.3 < v < 2.0 au. The choice is done to include the region of 
maximum effectiveness of the CTMC method: v > v e , with v e electron velocity. 

The presence of a second nucleus is clearly seen when one plots a versus i9 (Figure 3). One can 
see an increasing trend with i9, i.e. electron capture is favoured when the projectile impinges with 
a direction perpendicular to the molecule axis, plays here the same role of the angle (j> between 
the angular momentum of an aligned electron and the projectile direction in ion-atom collisions 
(see, for example, Figure 1 of Ref. [j/j]): with this parallelism in mind, the data may be compared 



with similar plots, for example, in |Mu22|. In comparison with those cases the effect is here much 
less marked, due to the fact that the electron probability distribution is smeared over a broader 
phase space volume. Nevertheless it seems possible to give at least a qualitative explanation of the 
trends in Figure 3 using propensity rules. As already explained in Section II one finds that, for a 
given value of the momentum p, the maximum of the probability of finding a matching between 
the velocities of the electron and the projectile is when p z = and-as p z = p cos ($)-this means 
1? = 7r/2 (see also Figure 1, where the electron distribution is localized close to p z = 0). cr(# = 90°) 
is almost constant up to v ~ 1 and only then falls down. In ref. jllj similar plots have been obtained 
for the scattering H + -H2 at high velocities ( > 1 MeV), within the Brinkman-Kramers formalism. 
There is discernible (see their Figure 6) a fluctuation, attributed to interference effects, which does 
not appear in our data (this was to be expected since, obviously, purely quantum mechanical effects 
cannot be included in our model). The same effect, even enhanced, is experimentally found in 
One way of looking at these data is plotting the anisotropy parameter 



versus v (Figure 4). a cx is <7; from Eq. (10) corresponding to the process of electron capture. A is 
oscillating but definitely assumes negative value, approaching zero while v increases. A < means 
that capture is favoured when the projectile and the molecular alignment are orthogonal. This is 
in agreement with other works (see, for example, the paper by Thomsen et al or that of Olson and 
Hoekstra in ref. where, furthermore, a more complex behaviour is also found, with changes of 
sign of A. 

In Figure 5(a) total cross sections (the same data of Figure 3, averaged over angle ■&) are shown as 
functions of impact velocity v. These data lend themselves to a comparison with structureless target 
scattering: in Janev [|3j it is empirically demonstrated how nucleus-hydrogen scattering follows a 
scaling law: the curve a cx /n 4 Z versus v 2 n 2 /Z ' 5 is universal, regardless of the initial principal 
quantum number n of the electron and of the charge Z of the nucleus. We may imagine to replace 
the diatomic molecule with a single particle, to which the electron is bound in a state defined by an 



effective (non integer) quantum number n e ff = 1/ . 2\E H +\ ~ 0.67 with our values. The agreement 



between our rescaled data and the universal curve by Janev yields an estimate of how much this 
modelling is justified. From Figure 5(b) one sees that the qualitative trend is the same, and the 
data are quite well interpolated by the fit in the middle of the range v. it is expected that, with 
increasing v, the electron-ion collisions closer and closer resemble two body processes, with a lesser 
influence of the target nucleus. In this situation the distribution function should not have influence. 
At the lower v's, the suggested fit underestimates the data; however, it is difficult to discerne how 
much of this discrepancy is due to the structure of the target and how much to the intrinsic defects 
of the CTMC method in this region of low energy. 

In order to have a further insight about the reliability of our results, we have compared them with 
previous calculations performed with other methods: in ref. ]l4| a calculation similar to ours has 
been carried on in the impact energy range from 100 keV to 5 MeV (2 < v < 14) using a distorted- 
wave model under different approximations: the simpler OBK approximation and the more refined 
correct-boundary-conditions Born serie (BIB) and the first order Bates series (Bal) (see and 
references therein for more details about these approximations). Figure 4 of ref. Jl4| shows the 
differential cross section for electron capture d<r/d(cosi?) as a function of i? at a collision energy of 
100 keV for iif + -_H^~ collisions. We have integrated the curves plotted and the results are shown 
in Figure 5(a). From this one may see that the accuracy of our calculation (at least for the single 
energy point available) is of the same order as the OBK approximation, and therefore overestimates 
the correct value, which should be close to that given by the BIB and Bal methods (which better 
fit the Janev' scaling law). 

Up to now, only g(p) has been taken into account to justify the results, so it is interesting 
to study the effects due to the spatial distribution g(r). Looking at the differential cross section 
da/db oc P(b)b for various impact energies and azimuthal angles, we noticed that the increase of a 
with & is due to the contribution from larger b's. This agrees with the results of jllj]. 

Finally, some words about the final state distribution. In our system about 80 — 90% of the total 
captures occur in the ground state. This is easily justified because the electron prefers to preserve 
its energy before and after the capture. A detailed study, looking for example at a dependence of 
this distribution from azimuthal angle or energy, would need a much larger amount of data, beyond 
the possibilities of the present study. 



A(v) 



<T cx {v,ti = 0°) +a cx (v,& = 90°) 



(11) 
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IV. SUMMARY AND CONCLUSIONS 



A series of numerical simulations has been performed on the charge-transfer collisions between 
protons and hydrogen molecular ions using classical methods. The interest of the subject relies on the 
comparison between this system and other, well studied, three-particle systems. Some conclusions 
which may be drawn from this study are: I) The CTMC method applied to this target is able to 
discerne its structure-as is seen from differential cross sections-but, with respect to quantal methods, 
its sensitivity is greatly reduced, as may be seen from the fact that no fluctuations due to interference 
effects are seen; II) Besides partial cross sections, also total cross sections seem to depend on the 
structure of the target, but this point is more difficult to stress since main differences appear at 
small v's, where CTMC is less reliable; III) The accuracy of the CTMC has been compared with 
quantal methods in the region of high v, limiting to total cross sections. It is found that -within the 
very small data set-the predictions of the CTMC well agree with those of the less refined versions 
of the quantum mechanical calculations, and slightly overestimate the more refined ones. 
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FIG. 2. Sketch of the geometrical arrangement of the collision. The two H+ are the nuclei of the molecule, p is the 
projectile and e~ the electron, zq is half the internuclear distance, 6 the impact parameter, $ the angle of impact and d the 
initial distance. For easiness, a coplanar collision is sketched: in the general case the projectile must be rotated out of the plane 
of an angle <j) ro t uniformly chosen within the range (0, 2ir). 
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FIG. 3. Electron capture cross section a cx versus azimuthal angle i? for different impact velocities. Errors bars are not 
shown as they are of the same size of the symbols. 
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FIG. 4. Anisotropy parameter A versus v. 
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FIG. 5. (a) Total electron capture cross section a cx , averaged over angle i), versus v. At v = 2 are also shown the data 
taken from ref. [14]: OBK approximation (full circle), BIB approximation (full square), Bal approximation (full triangle); see 



also the text, (b) The same data, but rescaled according Janev [23]: 
the present work, the solid line the fit from ref. [23]. 
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